DYNAMICS AND STATICS OF VORTICES ON 
A PLANE AND A SPHERE — I 



A. V. BORISOV 

Faculty of Mechanics and Mathematics, Department of Theoretical Mechanics 

Moscow State University 
Vorob'ievy gory, 119899 Moscow, Russia 
E-mail: borisov@uni.udm.ru 
A. E. PAVLOV 

Laboratory of Nonlinear Dynamics and Synergetics, Udmurt State University 
Universitetskaya Str. 1, 426034. Izhevsk, Russia 
E-mail: pavlov@uni.udm.ru 

Abstract 

In the present paper a description of a problem of point vortices on a plane 
and a sphere in the "internal" variables is discussed. The Hamiltonian equa- 
tions of motion of vortices on a plane are built on the Lie-Poisson algebras, 
and in the case of vortices on a sphere on the quadratic Jacobi algebras. The 
last ones are obtained by deformation of the corresponding linear algebras. 
Some partial solutions of the systems of three and four vortices are considered. 
Stationary and static vortex configurations are found. 

1 Dynamics of point vortices on a plane 
Introduction 

Let us consider the plane motion in boundless ideal liquid of the N rectilinear fila- 
ments with the strengths Tj, crossing the plane in the points with coordinates (xj, y^). 
It was shown by Kirchhoff [1] that the motion equations of such system can be writ- 
ten in the Hamiltonian form: 

TjXi = — , ViUi = -— , 1 ^ % ^ n, (1) 
% dxi 

with the Hamiltonian 

1 - 

# = -^E ,r * r > M ^ ( 2 ) 
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M i:j = (xi - xj) 2 + (y t - yjf, 
The Poisson brackets are of the form: 

N 



1 / 

i=i 1 ^ 



1 ( df dg df dg , 



dxi dyt dyi dx { 



The system of equations ((TJ) has, beside of energy (J2J), the first integrals due to 
the invariance of the Hamiltonian with respect to translations and rotations of the 
coordinate system: 



N N N 



i=l i=l i=l 

The set of integrals (j3J) is not in involution: 

A' 



{Q, P} = J2 r - i p > J > = W> J > = 2R ( 5 ) 



1=1 



Representation in "internal" variables 

In this paper another representation of problem is given. Let us choose the squares of 
mutual distances between the point vortices as new variables. For systems of N vor- 
tices the number of mutual distances is equal to N(N — 1)/2. Then, commuting My 
and Mm in the Poisson structure Q and using the Heron formulas, we as a result 
obtain some nonlinear Poisson brackets in the phase space determined by mutual 
distances. However, if we add new variables, the oriented areas of parallelograms 
spanned on triples of vortices with numbers i, j, k: 

Ay fc = (Tj - Ti) A (r k - r,i), r t = (x h y t ), (6) 

then in the phase space of the variables (My, Ay&, i, j, k — 1, . . . , N) the Lie- 
Poisson brackets appear. The dimension of the phase is made up of the number of 
mutual distances between vortices (the number of binomial combinations from N 
on 2) and the number of triangles formed by the triples of vortices (the number of 
binomial combinations of N things r at a time) (see Fig. 1): Cfj + Cf^ = C% +1 . 

Calculation of the Poisson brackets between these variables using (jHJ) leads to 
the following relations: 

{M^, M kl } = 4 (J^Sik - Y"$jk) A m + 4 - |T<^ Ay fe ; (7) 
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{Mij, A klm } = (JrSik - ^S jk j (M H - M im + M mj - M jt )+ 

+ (—5a - —Sji) (M mi - M lk + M kj - M jm )+ (8) 



(JrSim - -r<Wj (M kl - M u + M tj - M jk ); 



{Ajjfe, Ai mn } — —(Aj kn — Aj km ) + -—(Aj k i — Aj kn ) + — (Aj km — Aj kt )+ 

i i i i *• i 

+jr(^ikm — A ikn ) + -p— (A ikn — A ik i) + -^{Ai k i — A ikm )+ (9) 

-L — ( A —A N , 3 km / A _ A \,^kn(A A \ 

Thus, in the C^ +1 -dimensional phase space of the variables (M, A) the Lie- 
Poisson brackets are determined. The structure of this algebra is as follows: 

{M, M} c A, {M, A} c M, {A, A} c A. 

Besides the subalgebra A, this algebra contains subalgebras of problems of (N — 1), 
(JV- 2), . . ., 3 vortices. 

The Poisson brackets are degenerated. The linear Casimir function 

N 

A =Y1 r i r o M ij ( 10 ) 

is due to existence of the integral of angular momentum of vortex system. The 
linear Casimir functions correspond to the trivial geometric relations between the 
oriented areas of triangles. Let us write here, as an example, one such relation for 
a quadrangle Ay k + A ik i — Auj — Ay/. = (see Fig. 1). It is easy to see that the 
annihilator of the brackets is 

Fij k i = Ajjfc + A ik i — Auj — Aij k = 0. (11) 

The structure (JSJ) has also the additional Casimir functions. They came about 
from such relations as 

F ijk = {2A ijk f + M% + M% + Ml - 2(M ij M jk + M {j M ik + M jk M ik ) = (12) 

derived from the Heron formula, relating the areas and lengthes of sides of a triangle. 
In general even in the problem of four vortices the functions ()12j) are not Casimir 
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functions (they determine Casimir functions of subalgebra of three vortices). The 
detailed analysis and classification of Poisson brackets algebra Q we plan to present 
in the next part of our work (to be published). 

If we express from (|T2"|) and substitute in the equations of motion for the 
squares of mutual distances My, then we get the equations of Laura [2]. Having 
solved the equations of motion for the relative positions of vortices, we are able to 
find, using the quadratures and the initial conditions, the absolute coordinates of 
their positions on a plane at any time [3]. 



Partial cases of integrability of three and four vortices on a 
plane 

The problems of three and four vortices are of particular interest. The integrability 
of the problem of three vortices was indicated already by Poincare. The problem of 
four vortices is nonintegrable. 

Under reduction of the ten-dimensional algebra, corresponding to the problem 
of four vortices on the general level of six Casimir functions, the Poisson structure 
becomes nondegenerate. So, according to the Darboux theorem, there exist local 
canonical coordinates, in which the equations have the form of a usual Hamiltonian 
system with two degrees of freedom. For integrability of such system, besides the 
Hamiltonian, one more first integral is required, which, in general case, does not 
exist. 

Let us consider in detail a subalgebra ©-0, corresponding to the case of three 
vortices. In variables M, A it has the form: 

{M 1 , M 2 } = --^A 4 , {M 3 , M x } = -^A 4 , {M 2 , M 3 } = -^A 4 , 

1 3 1 2 1 1 

{M 1? A 4 } =(1--1^M 1 +Q- + ±-) (M 2 - M 3 ), 

1 1 \ (I 1 N (13) 



{M 2 , A 4 } = ( — - —j M 2 + ( — + — ) (M s - Mi,, 
{M 3 , A 4 } = ( i- - j-) Ms + Q- + ^ ) (M x - M 2) 



The rank of the algebra ()13j) is two. It corresponds on two-dimensional symplectic 
leaf to the system with one degree of freedom. Depending on the ratio of intensities, 
this algebra is a direct sum SO(3) ©for SO(l, 2) © R (the case Ti + T 2 + T 3 = 
is under consideration). Analysis of this problem will be also given in our next 
publication. 

One can point out at an interesting analogy between the problem of three vortices 
and the system of Lotka-Volterra, appearing in mathematical biology [4] . 
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The Hamiltonian of three vortices 



H = -^(T 2 T 3 lnM 1 + T 1 T 3 \nM 2 + T 1 T 2 lnM 3 ) (14) 
generates the phase flow: 

= ^ = s(sr si) (15) 

Expressing A4 from (|12|1 

A 4 = ±-^-Ml - Ml - Ml + 2(M l M 2 + MxM 3 + M 2 M 3 ), (16) 



we get the system: 



^ = ^mm; hh(hh - M ^ (17) 



If we introduce the regularizing time r: 

dr A, 



dt 2wM 1 M 2 M 3 
then we find that M, satisfy the Lotka-Volterra system 

d 



18) 



, M^IW^Ms-Ms). (19) 
dr 

Under transition of the system of vortices to a collinear configuration (A 4 = 0) 
a change of time flow is possible, so one must choose in (jTfij) the sign plus. Hence, 
we have demonstrated a piecewise trajectory isomorphism between two dynamical 
problems: the problem of three vortices and the three-dimensional Lotka-Volterra 
system, which takes place also in the problem of three vortices on a sphere, as it 
will be shown in the following section. 

Let us note a partial case of integrability of four- vortex problem [5]: 

4 

$^r i = o, Q = P = 0. (20) 

1=1 

In this case, the system has sufficient number of integrals Q, P, I, in involution (see 
the algebra of integrals (JSJ)). Its corresponding Casimir function is D on zero level: 
D = due to the relation between the first integrals [6]: 

D / 2 = fe^J J-Q 2 -P 2 . ( 21 ) 
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Additional invariant relations which are necessary for integrability in variables My 
one can find from the N(N — 1) identities [6] 



1 N 

fij =2^2 T ^{M jk - M ik ) = P{ Xl - Xj) + Q( Vi - Vj ), (22) 
fc=i 

being valid under ^2 i=1 1^ = 0. 

Let us indicate, for example, a particular solution of the system of four vortices 
if a central symmetry of the configuration takes place [3]. If Ti = r 3 , r 2 = and 
vortices form a parallelogram on a complex plane (— Z\ = z 3 , —z 2 = z±), then it 
preserves itself, i. e. in this case there are two invariant relations: 

M 1 = L x , M 3 = L 3 . (23) 

Having chosen an origin of coordinates, coinciding with the fixed symmetry center 
we can present the equations of motion in the form: 

dMi 



dr 
dM 3 

dr 



T 1 M 1 L 2 (M 3 - M 2 ) + T 2 M 1 M 2 (L 2 - M 3 ), 
T 1 M^L 2 {M 2 - Mi) + V 2 .\l 2 .\l,(.\l. - L 2 ), (24) 



^ = 2r 2 M 2 L 2 (M 1 - M 3 ), ^ = 2T 1 M 2 L 2 (M 3 - Mi). 
dr dr 

Here the notation: L 2 = M 24 is used and, similarily to f)18j) . the regularization 
of time r: 

£ = ^ (25) 

dt 2nM 1 M 2 M 3 L 2 v ; 

is made. The first pair of equations ()24j) describes motion of the sides of parallelo- 
gram and the second one describes motion of its diagonals. 
Equations (|2*4"|) are Hamiltonian with the Hamiltonian: 

H = 2rir 2 lnM 1 + r?lnM 2 + 2rir 2 lnM 3 + T 2 2 \nL 2 . (26) 

The Poisson brackets are quartic by elements of the algebra: 



{M 1 ,M 2 } = 


^M ± M 2 M 3 L 2 , 
1 1 


{M 1; M 3 } = 




l(Jr-fr) M X M 2 M 3 L 2 , 


{M 2 ,M 3 } = 


^M ± M 2 M 3 L 2 , 
1 1 


{M 1 ,L 2 } = 


~^M X M 2 M 3 L 2 , 
1 2 


{M 2 ,L 2 } = 


0, {M 3 ,L 2 } = ^M X M 2 M 3 L 2 
1 2 
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(27) 



The Poisson structure can be received from the bracket (|13|) under the restriction 
on the mentioned above invariant relations. Now, its rank is two. The geometric 
Casimir function is obtained from the relation between the squares of sides and the 
squares of diagonals in the parallelogram: 



For the real motions F is equal to zero. 

The Casimir function (fTU|) ( "angular momentum" of the system) takes the form: 



A reduction of system (|2~2j) on Casimir functions (J2~%)) . (J2T?)) reduces it to the 
one-dimensional integrable Hamiltonian one [6]. 

One must note that in contrast to the rigid body dynamics the equations of 
motion in the "internal" variables do not have a standard invariant measure as well 
as even a measure with the analytical positive density. As a result dynamics of these 
problems is deeply different. 

A reduction of system (j2U) on Casimir functions (J2~%)) . (ffilj) reduces it to the 
one-dimensional integrable Hamiltonian one [6]. 

One must note that in contrast to the rigid body dynamics the equations of 
motion in the "internal" variables do not have a standard invariant measure as well 
as even a measure with the analytical positive density. As a result dynamics of these 
problems is deeply different. 

2 Dynamics of point vortices on a sphere 
Introduction 

Let us consider a problem related to the previous one, a setting of which also goes 
back to the last century. In work [8] I. S. Gromeka obtained the equations of motion 
of point vortices on a sphere. The problem presents a simple model, describing 
a motion of vortex formations (cyclones etc.) in atmosphere of the Earth. The 
dynamical equations have been written in 1885 year yet were nearly forgot, though 
the problem in many respects is similar to the problem of motion of point vortices on 
a plane. In the papers [9, 10] there were pointed to integrability of three vortices on 
a sphere. In the papers [11, 12] there were proved the integrability of the restricted 
problem of four vortices. 

We shall show that the equations of motion of iV vortices on a sphere with 
the radius R can be presented in the Hamiltonian form with degenerated Poisson 
structure, assigned by Jacobi algebra JQ(C% +1 ) [13]. 




(28) 



D = V 1 M 2 + r 2 L 2 . 



(29) 
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In fact, in the usual canonical basis the equations of motion of iV vortices on a 
sphere can be written in the form: 

dH dH 
Qk = tt- > Pk = - , (30) 

opk oqk 

where the canonical coordinates q k , Pk are expressed through the spherical angles 9k, 
Pk and the strengths F k as 

q k = R\/fkip k , p k = R\/fkCOsd k . (31) 

The Hamiltonian function has the form: 

1 N 

H = -^J2 ' TiTk ln ( 4i?2 sin2 ^ fc / 2 )' ( 32 ) 

i,k=l 

where 

cos 7^ = cos 0i cos 6k + sin 9i sin 9 k cos(</?j — ifk) (33) 

is a cosine of the angle between i-th and fe-th radius vectors, connecting the center 
of sphere with the point vortices (see Fig. 2). 

The equations of motion of N point vortices in the spherical system of reference 
have the following form [9]: 



N 



0k 



1 ^ sin 0i sm((p k - (fi 



E, bill Ui blJ 
' ~cos lik 

1=1 



N 



. 1 ^ sin0 fc cos0i -cosflfcSinfliCOsOpfc f . 
Sm ^ = 4^ g ^ • (34) 



They are Hamiltonian equations with Poisson brackets 



{(Pi, cos6 k } = -^r- (35) 



The system (J3*UJ) always has four integrals of motion 

H = c , 



Fx = y ) Tj sin 9i cos (p { = ci, F 2 = 1^ sin ^ sin = c 2 , F 3 = Tj cos 0; = c 3 , 

iii 

(36) 

which, however, are not in involution. It is possible to demonstrate that the inte- 
grals F±, F2, F3 commute as components of angular momentum: 

{Fi, Fj} = EijkFk- (37) 
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Algebraic representation 

Let us introduce the "internal" variables: 

Mij = 4R 2 sin 2 7^/2, (38) 

being the squares of spans between arcs of corresponding vortices. The Hamiltonian 
depends only on introduced mutual variables: 



1 N 

# = -— ^T t r fe hiM 4fc . (39) 

i,k=l 

From the relations between the canonical coordinates qi, pj (J33j) it is possible to 
find commutators between the values M^: 

{M i:j , M kl } = 1 QrSik - ^6 jk j A ijt + i \J^5 U - i-tf^ A ijk , (40) 

where the following notations are used: 

A ijk = ri A Tj A r A . (41) 

Let us note that we use the same notations in both considered problems. Thus, the 
Poisson brackets between the values Mi k are proportional to external form of volume 
of parallelepiped spanned on the radius vectors of triple of vortices on a sphere. 
It is possible to demonstrate that the algebra of variables Mi k and Aij k is closed. 



{Mij, A klm } = R [ ^6 ik - ^5 jk ) (M u - M im + M mj - M jt )+ 



+R (Jrdu - ^-Sjij (M mi - M lk + M kj - M jm )+ 
+R [jrSim - jM;™) (M ki - M u + M u - M jk )+ 
+— \^5 ik - —S jk ] {M 3l M im - M mj Mu)+ 



2R \Ti r 7 



(42) 



+ jr \h 5il ~ r 6jl J ( M i mMik - M ^ M im)+ 
+ ^r (r 6im ~ r 6jm ' ( M i kMa ~ M ** M iO; 
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{Ajjfe, ^lmn\ 



8a 



+ ■ 



R(Aj kn Aj km 
R{^jkl — Aj kn 



+Y~ [ R( A,/,.,,, — A,/.., 



+ 



+ 



R\Aj, km Ai kn 
R(A ikn — A ik i 
R(A ik [ — A ikm 



4/ 
r 

Okm 



( R\^ijn Ajj 7; , 

J- fc 



-R(Ajj7 Ajj n 

A ij7 



+ p— ( R(A W ,, 



+ -^j^(M in Ajkm — M im Ajkn) ) + 



+ —(M u A jkn - M m A jM ) ) + 



+ —(M im A jkl - M u A jkm ) ) + 



+ -^^(Mj m Aikn — MjnAikm) ) + 



+ —(M jn A M - MjiA ikn ) ) + 



(43) 



+ -^j>(MjiAi km — Mj m A ik i) j + 
+ -^j^(M kn Aij m — M km Aij n )^j + 



+ —(M kl A ijn - M fcn A^)J + 
+ ■^j^(M km A i j l — M fc /Ajj m )^ . 



The quadratic algebra of vortices on a sphere is presented as a deformation of the 
linear algebra of vortices on a plane using the parameter — the radius of curvature of 
sphere R. The algebra of vortices on a plane (JTJ)-© and the corresponding Casimir 
functions of system on a sphere are obtained in the limit R — > oo : A v I R — > A, 
— > M from the algebra f)4U|) — . Here we symbolically marked by a symbol v 
ihe variables of the problem on the sphere. 

The linear Casimir functions of the quadratic algebra analogous to the corre- 
sponding functions of the linear algebra (fTU|) -()lH ) . The analog of the relations (fT^j) 
can be written as follows: 



Fijk =2A ijk 2 + R 2 {M 2 j + \ \% -r ivi jk) 



M 



-2R 2 (M i:j M jk + M ik M jk + MijM ik ) + M {j M jk M ki = 0. 



(44) 



Partial cases of integrabilily of three and four vortices on a 
sphere 

In the spherical problem, the analogue of the partial case on the plane is a situation 
when the integrals of motion (j3T?j) have zero values: 



0. 



(45) 
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The case is completely integrable because there are sufficient number of commuting, 
on zero level, integrals. 

It is not difficult to obtain identity, analogous to ({21} : 



N x 2 



D/2= UJ> -Fi-Fl-Fl (46) 



i=i 

Consequently, the considered case corresponds to the level of Casimir function: 

2 



D i 2 = ( R 52 T <) 



In order to find the invariant relations in variables it is possible to use identity, 
between absolute and mutual coordinates, analogous to 



1 N 

fij =2^2 T k( M jk - M lk ) = R{F!(xi - Xj ) + F 2 { yi - Vj ) + F 3 (zi - Zj )). (47) 



k=l 



3 Stationary vortex states 

The obtained forms of the dynamical equations on a plane and on a sphere can 
be successfully used for finding the steady-state configurations, being the partial 
solutions of the equations of motion. 

Requirements of N(N — l)/2 mutual distances My = being constant are the 
conditions of stationary state of vortex configurations. They have the same form 
both for the configurations on a plane and on a sphere in the "internal" variables as 
its Hamiltonians (J2J), an d commutation relations {M^, M k {\ (J7|), have the 
same form: 

N f 1 1 \ 

(4£ 



ir'frr- tf-) r/A«,(M) = 0. 



The conditions ()48j) are very suitable for finding symmetric configurations. On 
the plane there is solution of N vortices of equal strengths T, disposed in vertices of 
regular polygon, refined into a circle of radius Rq. The system rotates with angular 
velocity 

n-W- 1 *. (49) 

The analysis of stability of such configurations in linear approximation was carried 
out in [15, 16]. 
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For investigation of such configurations J. J. Thomson was awarded with Adams prize 
in 1883 [15]. These solutions were at the basis of the theory of vortex atoms pro- 
pagandized by W.Kelvin (before quantum mechanics). Analogous configuration of 
identical vortices on a sphere is disposed, according to (J52j) . on latitude 6 = 9 , 



The angular velocity of the chain of vortices decreases from poles to equator and 
on the equator the configuration is static. 

The next solution of the system (|48jl are collinear configurations: N identical 
vortices on a line rotate with some angular velocity around axis perpendicular to 
the plane. On the sphere of radius R the collinear configuration ( N atoms on the 
meridian) rotates around the axis z with angular velocity Q (see Fig. 4). 

Collinear configurations of vortices on a sphere were not discussed in literature 
before. We can find them using the following conditions: coordinates (fi of the point 
vortices have the equal values or differ on tt. In addition: Aj^(M) =0. It is easy 
to notice that coordinates 6^, (0 ^ 6 ^ 2n) are roots of the following system of 
trigonometric equations: 



i=l v 7 

The chain of vortices (JSTj) deserves more attention. The roots of the system 
assigned the collinear configuration 9i can be found as equilibrium position of the 
system of N particles on a circle set by the Hamiltonian 



Evidently, the positions of equilibrium of the chain of particles coincide with roots 
of the system Such connection in plane case was marked by Calogero [17]. 

The locations of identical vortices on a line are set by zeros of Hermitean polyno- 
mials of N-th power. The system (J52j) under Q = was considered in [18], where 
statistical properties of levels of energies of one-dimensional classical Coulomb gas 
are considered. In the paper [19] the results of equilibrium analysis of the system 
(J52J) under condition ft — : 0j = 9q + kit j N, k = 1, ... , N are adduced. 
The problem defined by the Hamiltonian (J52j) is usually nonintegrable. 
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Figure 4: 
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